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PROCEDE DE C&R&CTERI SATICM D'XM SIGML SOKTORE 

L 1 invention concerne un procede de caracterisation 
selon des parametres specif iques, d'un signal sonore 
evoluant selon le temps dans differentes bandes de 
frequences . 

5 Le dornaine de 1 ' invention est celui de la 

reconnaissance de signaux sonores appliqu^e en 
particulier a 1 1 identification d'oeuvres musicales 
utilis<§es sans autorisation. 

En effet, le developpement des techniques de 

10 nuni^risation et du multimedia a provoque une 
augmentation considerable de telles utilisations 
frauduleuses. II en resulte une nouvelle difficulte 
pour les organismes charges de collecter les droits 
d'auteur, puisqu'il faut etre en mesure d' identifier 

15 ces utilisations en particulier sur les reseaux 
num^riques interactifs tel Internet, afin d'effectuer 
de fagon satisf aisante la perception et ,1a repartition 
des remunerations dues aux auteurs de ces ceuvres 
musicales . 

20 Dans la suite, pour ne pas se limiter aux oeuvres 

musicales, on considerera plus generalernent un signal 
sonore. 

Le but de la pr^sente invention est done de 
constituer une base de donnees de signaux sonores, 
25 chaque signal sonore 6tant caract6ris£ par une 
empreinte de sorte que, etant donn6 un signal sonore 
inconnu que 1 1 on caracterise de la meme fagon, on 
puisse effectuer une recherche et une comparaison 
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rapides de l'empreinte de ce signal inconnu avec 
1' ensemble des empreintes de la base de donnees . 

L'empreinte est constitute de parametres 
specif iques determines de la fagon suivante. 

On decompose dans un premier temps le signal 
sonore dont 1 ' amplitude x(t) varie avec le temps t, 
selon differentes bandes de frequences k : x(k,t) est 
1' amplitude du signal sonore f iltre - dans la bande de 
frequence k et represents figure la) . 

On calcule, comme represents figure lc) 1 ' energie a 
court terme E(k,t) de ce signal sonore f iltre a l'aide 
d'une fenetre h(t) representee figure lb), ayant un 
support de 2N secondes . On reitere ce calcul en faisant 
glisser cette fenetre toutes les S secondes. 

Ces valeurs E(k,t) constituent des parametres 
specif iques d'un extrait de 2N secondes du signal 
sonore x(k,t) dans la bande de frequences k. 

On peut obtenir d'autres parametres en calculant 
pour differentes bandes de frequences j, 1' energie de 
E(k,t) a l'aide d'une fenetre h' (t) representee figure 
2b) , ayant un support de 2N' secondes ; on reitere ce 
calcul en faisant glisser cette fenetre toutes les S' 
secondes : on obtient F(j,k,t) represents figure 2c). 
On normalise les valeurs F(j,k,t) par rapport a leur 
maximum de fagon a les rendre independantes de 
1' amplitude du signal sonore. 

Ainsi normalisees, ces valeurs constituent des 
parametres specif iques d'un extrait de 2N' secondes du 
signal sonore x(k,t) dans la bande de frequences k. 

On peut egalement calculer la phase de E(k,t) pour 
differentes bandes de frequences j : on obtient 
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P(j,k,t). On normalise les valeurs P(j,k,t) par rapport 
k une valeur de reference P(l,j,t) et on obtient alors 
d'autres parametres specif iques d'un extrait de 2N' 
secondes du signal sonore. 
5 D'autres parametres peuvent §tre ajout^s tels que 

la valeur moyenne de 1 * energie E(k,t). 



L' invention a pour objet un procede de 
caracterisation selon des parametres specif iques, d'un 

10 signal sonore x(t) evoluant selon le temps t pendant 
une duree D dans differentes bandes de frequences k et 
alors note x(k,t), principalement caracteris^ en ce 
qu'il consiste a m^moriser le signal x(t) , a calculer 
1' energie E(k,t) dudit signal x(k,t) pour chacune 

15 desdites bandes de frequences k, k variant de 1 a K et 
selon une fen§tre temporelle h(t) d'une duree 2N, 
meraoriser les valeurs de 1 1 energie E(k,t) obtenues, ces 
valeurs constituant les parametres specif iques d'un 
extrait d'une dur^e 2N du signal sonore x(t) et k 

20 reiterer ce calcul a intervalles reguliers pour obtenir 
1' ensemble des parametres specif iques pour la duree D 
du signal sonore x(t) . 

II consiste en outre & calculer et memoriser 
1' Energie F(k,j,t) de E(k,t) pour des bandes de 

25 frequences j , j variant de 1 a J, selon une fenetre 
temporelle h' (t) d'une dur6e 2N' , les JxK valeurs de 
1 ' Energie F ( j , k, t) obtenues constituant les parametres 
specif iques d'un extrait d'une duree 2N' du signal 
sonore x(t) et a reiterer ce calcul a intervalles 
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reguliers pour obtenir 1' ensemble des parametres 
specif iques pour la duree D du signal sonore x(t) . 

II consiste eventuellement a calculer la phase 
p(j,k,t) de l'energie E(3c,t) pour des bandes de 
frequences j, j variant de 1 a J, j different de k, et 
a inclure les valeurs de la phase P(j,k,t) obtenues 
parmi les parametres specif iques du signal sonore x(t) . 

II peut aussi consister a calculer pour chaque 
bande de frequence j, la valeur moyenne de l'energie 
E(k,t) sur 2N' secondes, a reiterer ce calcul a 
intervalles reguliers pour obtenir 1 ' ensemble des 
parametres specifiques pour la duree D du signal sonore 
x(t) et k inclure les valeurs moyennes obtenues parmi 
les parametres specifiques du signal sonore x(t) . 

Selon une caracteristique, il consiste k considerer 
les parametres specifiques d'un signal sonore x(t) 
comme les composantes d'un vecteur representatif de 
x(t) , a positionner les vecteurs dans un espace a 
autant de dimensions que de parametres, a definir des 
classes regroupant les vecteurs les plus proches et a 
enregistrer lesdites classes . 

Les classes presentant des distances inter-classes 
et des distances intra-classes, le procede consiste 
avantageusement a selectionner parmi les parametres 
specifiques, les parametres permettant d' obtenir des 
distances inter-classes relativement grandes devant les 
distances intra-classes et a enregistrer les parametres 
selectionnes. 

L' invention concerne egalement un dispositif 
d' identification d'un signal sonore, caracterise en ce 
qu'il comprend un serveur de base de donnees comprenant 
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des moyens pour la mise en oeuvre du procede de 
caract^risation d'un signal sonore selon des parametres 
specif iques, tel que decrit pr^cedemment et des moyens 
de recherche dudit signal sonore dans la base de 
5 donnees . 

De preference, les moyens de recherche comprennent 
des moyens de reconnaissance directe de la classe k 
laquelle ledit signal sonore appartijent et des moyens 
de recherche de la classe par comparaison des 
10 parametres specifiques du signal sonore inconnu avec 
ceux de la base de donnee, la classe etant choisie, par 
exemple, par la m^thode des plus proches voisins. 



15 D'autres particularity et avantages de 1 ' invention > 

apparaitront clairement k la lecture de la description }k 
faite a titre d 1 exemple non limitatif et en regard des ; 
dessins annexes sur lesquels : 

les figures la), lb) et lc) representent 

20 respect ivement des courbes schematiques de variation 
d'un signal sonore x(ki,t) filtre dans une bande de 
frequences k i7 d'une fenStre de Hamming h(t) et de 
1'energie E(ki,t) a court terme du signal x(ki,t), 

les figures 2a), 2b) et 2c) representent 

25 respectivement des courbes schematiques de variation de 
1'energie E(ki,t) pour la bande de frequence k±, d'une 
fen§tre de Hamming h' (t) et de l'energie F(j«,ki,t) de 
E(k ir t) pour la bande de frequences j m , 

la figure 3 illustre schematiquement un ensemble de 

30 vecteurs V[x(t)] constituant l'empreinte d'un signal 
x(k, t) , 
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la figure 4 illustre schematiquement le stockage 
d ' empr eintes , 

la figure 5 represente une classification des 
signaux sonores selon deux parametres, 

la figure 6 illustre une methode de recherche d'un 
signal sonore par la methode du plus proche voisin. 

la figure 7 represente sch&natiquement un serveur 
de base de donnees de stockage des empreintes de 
signaux sonores . 



Les signaux sonores que l'on traite selon ce 
procede de caracterisation sont des signaux sonores 
enregistr^s, notamment sur des disques compacts 
(Compact Disk en anglais) . 
15 on va considerer par la suite que le signal sonore 

x(t) est un signal numerique <§chantillonn6 & une 
frequence d' echantillonnage fe, par exempl'e 11025 Hz 
correspondant au quart de la frequence 
d' echantillonnage courante des disques compacts qui est 

20 de 44100 Hz. 

On peut cependant caracteriser un signal sonore 
analogique : il faut au prealable le convertir en un 
signal numerique au moyen d'un convertisseur 
analogique-num^rique . 

25 Le signal sonore x(k,t) represente figure la) pour 

k=ki est done un signal numerique echantillonne a la 
frequence fe et obtenu apres un filtrage dans une bande 
de frequences ki. Chaque valeur de ce signal numerique 
echantillonne est cod6e par exemple sur 16 bits. Les 

30 bandes de frequences sont des bandes du spectre audible 
variant d' environ 20 Hz a 20kHz et decoupe en K {k 
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varie de 1 a K) bandes de frequences, K=127 par 
exemple . 

L'energie a court terme E(k, t) representee figure 
lc) pour k=ki, est calcul<§e sur une fen@tre h(t) de 2N 

5 secondes, par exemple une fenetre de Hamming ayant un 
support d 1 environ 23 ms representee figure lb). 

E(k,t) est le carre du module d'une transformee du 
signal sonore echantillonne x(t) dans le plan temps- 
frequence ou dans le plan temps-echelle. Parmi les 

10 transformees utilisables, figurent la transformee de 
Fourrier, la transformee en cosinus, la transformee de 
Hartley et la transformee en ondelettes . Un banc de 
filtres passe-bandes realise egalement une telle 
transformee. La transformee de Fourier a court terme 

15 permet une representation temps -frequence adaptee k 
1 'analyse des signaux musicaux. Ainsi l'energie E(k,t) 
s ' 6crit : 

2 



E(k,t) = 



2 x(t + n/fe).h(n/fe).e 

n=-N 



20 avec i tel que i 2 = -1 

On fait glisser la fenetre sur le signal sonore 
toutes les S secondes, par exemple toutes les 10 ms. 
E(k,t) sera ainsi echantillonne toutes les 10 ms : on 
aura E(k,t 0 ) puis E(k,ti) avec ti=t 0 +10 ms etc. 

25 On code ainsi toutes les S secondes le signal 

sonore x(t) par un vecteur a K composantes E(k, t) , 
chacune de ces composantes codant l'energie de 23 ms du 
signal sonore x(t) dans K bandes de frequences. 
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On obtient d'autres parametres en reproduisant en 
quelque sorte les calculs precedents et en les 
appliquant cette fois a E(k,t) comme represente figure 
2a) a figure 2c) . 

5 On filtre 1 ' energie E(k,t) dans J differentes 

bandes de frequence : E(j,k,t) est 1' energie. E(k,t) 
f iltree dans la bande de frequences j , j variant de 1 a 
J avec par exemple J=51 . 

On calcule alors F(j,k,t) representee figure 2c) 

10 pour k=ki et j=j m , sur une fenetre b' (t) de 2N' 
secondes, par exemple une fenetre de Hamming ayant un 
support de 10 s. Ainsi, on peut ecrire, avec i tel que 
i 2 = -1 : 

' 2 EOC t + n/fcO-h'Cn/feO.e- 4 ^' 



15 



F(j,k,t) = 



n=-N" 



Dans notre exemple, toutes les secondes (S'=l), le 
signal sonore x(t) est code par 127x51 parametres 
F{j,k,t), chaque reel F(j,k,t) representant 1 ' energie 
de dix secondes (2N'=10) du signal d" energie E(k, t) 

20 dans la bande de frequence j . 

De maniere a rendre F(j,k,t) independant de 
1' amplitude du signal qui peut etre plus ou mo ins fort, 
on rapporte ces valeurs a une valeur de reference, en 
1' occurrence la valeur maximale de F M (j,k,t) pour tous 

25 les k et j consideres. On obtient ainsi KxJ parametres 
F(j,k,t) /F M (j,k,t) . 
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On calcule en outre toutes les 2N' secondes, la 
phase de l'energie E(k,t) dans chacune des bandes de 
frequences j: P(j,k,t). 

Pour cela, on calcule 1 'argument de la transformee 
de Fourrier de E(k,t) dans chacune des bandes de 
frequence j : 



Comme precedemment, on rapporte ces valeurs a une 
valeur de reference, en 1" occurrence la valeur de 
P(j / k / t) pour la seconde bande de frequences (j=l) 
consider6e car la reference temporelle de 1 'echantillon 
est inconnue : on ne'connait pas l'origine des temps. 
Pour cela, on calcule les phases rapportees (p(j,k,t) 
par les formules suivantes : 
9(l,k,t) = P(l,k f t) 

<p(j,k,t) = P(j,k,t) - P(l,k,t). f (k) , pour k>l. 



ou les f (k) sont les frequences centrales des canaux k. 
On obtient ainsi KxJ parametres correspondant aux 
valeurs de la phase rapportee <p (j,k,t). 

On peut ^galement considerer d'autres parametres 
notamment, les valeurs moyennes de l'energie E(k,t) sur 
2N' secondes et ce pour chaque bande de frequences j : 
E(j,k,t) . 




f (1) 
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L' ensemble de ces parametres normalises 
definissent a intervalles reguliers, une empreinte que 
l'on peut considerer comme un vecteur V(x(t)). 
L 1 ensemble des parametres normalises par exemple 
5 F<j,k,t)/F M et P{j,k,t)-P(j,l,t) definissent toutes les 
S' secondes, une empreinte que l'on peut considerer 
comme un vecteur V(x(t)) a 2xKxJ dimensions (2x127x51 
soit environ 13 000 dans notre exemple), une dimension 
par parametre, chaque vecteur caracterisant un extrait 
10 de 2N' secondes du signal sonore x(t), 10 secondes dans 
notre exemple . 

On reitere cette caracterisation toutes les S' 
secondes, toutes les secondes par exemple (S'=l). 

Comme represents figure 3, un signal x(t) durant T 
15 secondes est finalement caracterise par L vecteurs V, L 
etant environ egal a T/S' . 

Pour un signal sonore durant 10 mn soit 600s, on 
obtient 600 vecteurs c'est-^-dire 600x2x«JxK parametres. 

Ces vecteurs sont memorises dans la zone de 
20 stockage 10 d'une base de donnees hebergee sur un 
serveur ou sur un disque compact. On a reprSsente 
figure 4 1 ' ensemble des vecteurs V d'un signal ou d'une 
oeuvre A par VA, de meme VB pour une oeuvre B, etc. 

II est souhaitable de reduire le nombre de 
25 composantes de ces vecteurs c'est-a-dire le nombre de 
parametres de facon a obtenir un vecteur ou une 
empreinte de taille plus reduite en vue de son stockage 
dans la base de donnees. De plus, lorscru'il s'agira de 
comparer 1 ' empreinte d'un signal sonore inconnu avec 
30 celles de la base de donnees, il sera souhaitable que 
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le nombre de param^tres a comparer soit reduit pour que 
cette recherche soit effectuee rapidement. 

Or ces param£tres n'apportent pas tous la m§me 
quantity d' information, certaines pouvant etre 
5 redondantes voire inutiles. C'est pourquoi on 
selectionne les parametres les plus pertinents parmi 
tous les pajrametres, par un calcul d' information 
mutuelle presente notamment dans la . publication PROC. 
ICASSP'99, Phoenix, Arizona, USA, Mars 1999 H.YANG, 
10 S.VAN VUUREN, H.HERMANSKY, « Relevancy of Time- 
Frequency Features for Phonetic Classification Measured 
by Mutual Information ». On limite ainsi K a Ki et J & 
Oi. 

On va presenter une methode de selection de ces 
15 parametres . r : , 
Chacune des empreintes de ces signaux sonores V 
c'est-&-dire chacun de ces vecteurs est classe dans un i 
espace R a N dimensions, N etant le nombre de * 
composantes des vecteurs. Pour des raisons de 
20 simplification, un exemple de classification pour des 
vecteurs a 2 dimensions PI et P2 est represents figure 
5. 

On def init des classes C (ra) regroupant les 
vecteurs par proximite, m variant de 1 a M. On peut par 
25 exemple decider qu'une classe correspond a une ceuvre 
musicale : dans ce cas M est le nombre d'oeuvres 
musicales memorisees dans la base de donnSes, 

II resulte du calcul de 1 ' information mutuelle 
entre ces classes C (m) et les parametres, que la 
30 pertinence des parametres est liee aux distances inter 
et intra classes : des parametres pertinents 
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garantissent des distances inter-classes d relativement 
grandes comparees aux distances intra-classes D. 

En ne retenant que les param£tres pertinents, on 
definit ainsi Ki et Ji. 

On peut par exemple considerer cinq (K x =5) bandes 
de frequences respectivement centrees sur 344 Hz, 430 
Hz, 516 Hz, 608 Hs et 689 Hz. 

Des essais ont ete effectuts en prenant Ji = 3. 

Les classes C (m) sont alors constitutes a partir 
des vecteurs V q (x) ne comportant plus que 2xKi.xJi 
composantes . 

On va donner un exemple pour K x =5 et Ji=3, de la 
taille memoire d'une base de donntes contenant 1000 
heures de musique et en considerant comme parametres 
E(k,t) et F(j,k,t), chacun de ces parametres etant code 
sur 4 octets. 

Les parametres E(k,t) .calcules toutes les 10 ms 
occupent 1000x3600x100x5x4 octets soit environ 7 Giga 
octets . 

Les parametres F<j,k,t) calcules toutes les 
secondes occupent 1000x3 600x3x5x4 octets soit environ 
200 Mega octets. 

Ces parametres sont associes aux references des 
signaux sonores ; si 1 1 on considere que les references 
contiennent 100 caracteres codes chacun sur un octet, 
ces references occupent 1000x10x100 octets soit environ 
1 Mega octet. 

Une telle base de donnees occupera finalement 
environ 7 Giga octets. 
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Lorsqu'on souhaite identifier un signal sonore 
inconnu, on en etablit d'abord 1' empreinte, ref^rencee 
V(xinc) a la figure 6, comme d^crit pr^cedemment , 
sachant que le signal sonore inconnu peut dtre une 
5 oeuvre musicale complete ou un extrait de cette oeuvre. 

La recherche de la classe de cette empreinte dans 
la base de donnees consiste alors selon une methode 
classique illustree figure 6, a comparer les parametres 
de cette empreinte V(xinc) k ceux des empreintes de la 
10 base de donnee. Les empreintes les plus proches, dites 
les plus proches voisins, definissent la classe de la 
fagon suivante : la classe est celle de la majorite des 
plus proches voisins. 

15 Un serveur de base de donnees 1 est 

schematiquement represents figure 7 . II comprend une 
zone de stockage 10 des donnees de la* base dans 
laquelle sont memorisees les empreintes des signaux 
sonores assorties de leurs references, II comprend en 

20 outre une m&noire 11 dans laquelle sont stock^s les 
programmes de caracterisation et de recherche 
precedemment decrits, un processeur 12 assorti de 
memoir es de travail pour mettre en oeuvre ces 
programmes. II comprend bien sftr une interface 

25 d' entree-sortie 13 et un bus 14 reliant ces divers 
elements entre eux. 

Lorsqu'il s 1 agit de rentrer dans la base de 
donnees 1, de nouveaux signaux sonores, 1 1 interface 13 
regoit le signal x(t) assorti de ses references ; s'il 

30 s'agit seulement d'un signal inconnu & identifier, 
1" interface 13 regoit seulement le signal x{t) inconnu. 
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En sortie, 1 ' interface 13 fournit vine reponse a la 
recherche d'un signal inconnu. Cette reponse est 
negative si le signal inconnu n'existe pas dans la zone 
de stockage 10 ; si le signal a ete identifie, la 
reponse comporte les references du signal identifie. 
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1 . Procede de caracterisation selon des parametres 
specif iques, d ! un signal sonore x{t) evoluant selon le 
temps t pendant une duree D dans differentes bandes de 
frequences k et alors not£ x(k,t), characterise en ce 

5 qu'il consiste k memoriser le signal x(t) , a calculer 
1 1 energie E(k,t) dudit signal x(k,t) pour chacune 
desdites bandes de frequences k, k variant de 1 a K et 
selon une fen§tre temporelle ti(t) d'une duree 2N, 
memoriser les valeurs de 1* energie E(k,t) obtenues, ces 
10 valeurs constituant les parametres specif iques d'un 
extrait d'une duree 2N du signal sonore x(t) et a 
reiterer ce calcul k intervalles reguliers pour obtenir 
1' ensemble des parametres specif iques pour . la duree D 
du signal sonore x(t) . 

15 

2. Procede selon la revendication prec^dente, 
caracterise en ce qu'il consiste a calculer et 
memoriser 1 ' energie F(k,j,t) de E(k,t) pour des bandes 
de frequences j , j variant de 1 a J, selon une fenetre 

20 temporelle h' (t) d'une dur6e 2N' , les JxK valeurs de 
l'^nergie F(j,k,t) obtenues constituant les parametres 
specif iques d'un extrait d'une duree 2N 1 du signal 
sonore x(t) et a reiterer ce calcul a intervalles 
reguliers pour obtenir 1' ensemble des parametres 

25 specif iques pour la duree D du signal sonore x(t) . 
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3. Procede selon l'une des revendications 
precedentes, caracterise en ce qu'il consiste a 
calculer la phase rapportee (p (j , k f t) de 1 1 energie 
E(k,t) pour des bandes de frequences j, j variant de 1 

5 a J, et k inclure les valeurs de la phase rapportee 
<p(j,k,t) obtenues parmi les parametres sp6cif iques du 
signal sonore x(t) . 

4. Procede selon l'une des revendications 
10 precedentes, caracterise en ce qu'il consiste a 

calculer pour chaque bande de frequence j , la valeur 
moyenne de 1 1 energie E(k,t) sur 2N'' secondes, k 
reit^rer ce calcul a intervalles reguliers pour obtenir 
1' ensemble des parametres specif iques pour la dur£e D 
15 du signal sonore x(t) et k inclure les valeurs moyennes 
obtenues parmi les parametres specif iques du signal 
sonore x(t) . 

5. Procede selon l'une des revendications 
20 precedentes, caracterise en ce qu'il consiste k 

considerer les parametres specif iques d'un signal 
sonore x(t) corame les composantes d'un vecteur 
representatif de x(t), k positionner les vecteurs dans 
un espace k autant de dimensions que de parametres , a 
25 definir des classes regroupant les vecteurs les plus 
proches et k enregistrer lesdites classes. 

6. Procede selon la revendication precedente, 
caracterise en ce que les classes presentent des 

30 distances inter-classes et des distances intra-classes 
et en ce qu'il consiste a selectionner parmi les 
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param£tres specif iques, les parametres permettant 
d'obtenir des distances inter-classes relativement 
grandes devant les distances intra-classes et a 
enregistrer les parametres selectionn^s . 

5 

7, Dispositif d' identification d'un signal sonore, 
caracterise en ce qu'il comprend un serveur.de base de 
donnees comprenant des moyens pour la mise en oeuvre du 
procede de caracterisation d'un signal sonore selon des 

10 parametres specif iques, selon l'une quelconque des 
revendications prec^dentes et des moyens de recherche 
dudit signal sonore dans la base de donnees* 

8. Dispositif selon la revendication precedente 
15 prise en combinaison avec la revendication 5 ou 6, 

caracterise en ce que les moyens de recherche 
comprennent des moyens de reconnaissance de la classe a 
laquelle ledit signal sonore appartient et des moyens 
de comparaison, par la methode du plus proche voisin, 
20 des parametres specif iques du signal . sonore inconnu 
avec les parametres specif iques de la base de donnees • 
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